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We derive the class of covariant measurements which are optimal according to the maximum 
likelihood criterion. The optimization problem is fully resolved in the case of pure input states, under 
the physically meaningful hypotheses of unimodularity of the covariance group and measurability 
of the stability subgroup. The general result is applied to the case of covariant state estimation for 
finite dimension, and to the Weyl-Heisenberg displacement estimation in infinite dimension. We also 
consider estimation with multiple copies, and compare collective measurements on identical copies 
with the scheme of independent measurements on each copy. A "continuous- variables" analogue of 
the measurement of direction of the angular momentum with two anti-parallel spins by Gisin and 
Popescu is given. 
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I. INTRODUCTION 



State estimation is a unique kind of quantum measure- 
ment in the quality of information that it provides. In 
fact, the knowledge of the state of a quantum system 
enables the evaluation of any ensemble average, which is 
equivalent to the possibility of performing any desired ex- 
periment on the system. For its intrinsic versatility such 
unconventional type of quantum measurement is of in- 
terest for the new technology of quantum information ^ 
in the estimation of parameters that do not correspond 
to observables — such as the phase of an electromag- 
netic field |3] — but also as a method to achieve quantum 
cloning |Q, 5|, whence in designing eavesdropping strate- 
gies for quantum cryptography 

An exact state estimation without any prior knowledge 
on the form of the state is impossible due to the no- 
cloning theorem |^. This also reflects the fact that 
an optimal approximate state estimation would not be 
achievable as an orthogonal measurement, since the state 
estimation is a kind of " informationally complete" mea- 
surement Ip^ . More generally, one can have some prior 
knowledge of the form of the state, i.e. by parameterizing 
it with a restricted set of variables. This is the typical sit- 
uation of the Quantum Estimation Theory of Helstrom 
where the goal is to determine a multidimensional 
parameter of a state transformation. When the set of 
states to be discriminated are orthonormal the parame- 
ter corresponds to an " observable" whose eigenstates are 
the set itself, and the estimation is exact. However, in 
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practice it happens very often that the multidimensional 
parameter cannot be described by an observable (e.g. it 
is a phase of a field, or it corresponds to a set of non com- 
patible observables), whence a measurement represented 
by a so-called positive operator valued measure (POVM) 
needs to be performed. 

For a state estimation which is not equivalent to the 
measurement of an observable we have a choice of in- 
finitely many POVM's achieving the same task with dif- 
ferent strategies. Indeed there is no universal criterion 
which is optimal for all situations, and one needs to de- 
fine the appropriate figure of merit pertaining to the par- 
ticular problem. Once the optimization problem is solved 
in terms of an optimal POVM, one can then address the 
problem of the feasibility of the measurement appara- 
tus by classification of orthogonal dilations of the POVM 
, or else compare the performance of actual de- 
vices to the ultimate theoretical limit. 

A statistically meaningful optimization strategy is the 
maximization of the likelihood that the true value of the 
estimated parameter coincides with the outcome of the 
measurement. Such a strategy is actually very general, 
since for measurements which are group- covariant, op- 
timization of a generic goal function corresponds to op- 
timization of the likelihood for a different input state. 
Physically, " group-covariance" means that there is a 
group of transformations on the probability space which 
maps events into events, in such a way that when the 
quantum system is transformed according to one element 
of the group, the probability of the given event becomes 
the probability of the transformed event. This situation 
is very natural, and occurs in most practical applications. 
For example, the heterodyne measurement Q is co- 
variant under the group of displacements of the com- 
plex field, which means that if we displace the state of 
radiation by an additional complex averaged field, then 
the output photo-current will be displaced by the same 
complex quantity. Other examples of covariant measure- 
ments are the quantum estimation of a " spin orientation" 
p5| , p^ , pT| |l8| , or of the phase shift of an electromagnetic 
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field |, |T|,|20|. 

The statistics of the measurement can be improved by 
using many copies of the same quantum system. In this 
scenario, it is relevant for experiments to distinguish the 
measurements achievable by local operations and classi- 
cal communication (LOCC) from more general schemes 
that require entanglement. Unfortunately, a useful clas- 
sification of LOCC schemes is still missing. Alterna- 
tively, one can give just a mathematical categorization 
in terms of the POVM of the measurement: i) "indepen- 
dent" measurements, corresponding to tensor product 
of independent POVM's; ii) "separable" measurements, 
corresponding to POVM's where each element is separa- 
ble; iii) "non separable" or "entangled" measurements, 
corresponding to POVM's where some element is entan- 
gled. In the first category measurements are performed 
independently on each copy. In the separable class, on 
the other hand, the measurement can be performed by 
means of separable operations, hence all LOCC schemes 
are included in this category. Notice, however, that not 
all separable operations can be implemented locally (see, 
e.g. , the case of nonlocality without entanglement of Ref. 

Finally, the class of entangled POVM's represents 
the most general scheme of measurement, and opens the 
exponential growth of the Hilbert space dimension versus 
the number of copies N, with the possibility of largely 
surpassing the statistical efficiency of the independent 
measurement schemes |3| ^ . However, as already 

noticed in Ref. [ psf , for the maximum likelihood strategy 
the optimal schemes can be surprisingly achieved by sep- 
arable measurements, and here we address this issue for 
covariant measurements. Under the general assumption 
of square-summable representation we derive a general 
"canonical form" for the optimal measurements for pure 
input states, corresponding to a POVM which is separa- 
ble or entangled, depending on the group representation. 

After introducing in Section || the precise formulation 
of the covariant state estimation problem, in Section III 



we derive some useful mathematical identities for group 
integrals which are then used to algebraically character- 
ize covariant measurements. This also helps us in deriv- 
ing a simple upper bound for the maximum likelihood in 
Section [V , along with the canonical form of the optimal 
measurement given in terms of the group representation. 
Examples of the canonical form are given in Section m 
dimension < cxd for the group SU{d) — corresponding 
to the estimation of an unknown pure state — and in infi- 
nite dimensions for the estimation of displacements on 
the phase space. The case of multiple copies is then 
analyzed, discussing the occurrence of entangled versus 
separable POVM's. For the estimation of displacements 
on the phase space, the case of two copies experienc- 
ing opposite shifts in momentum is also analyzed — the 
continuous- variables analogue of the measurement of di- 
rection of the angular momentum with two antiparallel 
spins by Gisin and Popescu For coherent states it 
is shown that such a scheme provides a better estimation 
of the displacement as compared to the conventional case 



of identical displacements. 



II. THE PROBLEM 

Whenever a quantum system S undergoes a physical 
transformation belonging to a group G, its state is trans- 
formed according to an appropriate representation of G 
on the Hilbert space H of the system S. In the following, 
we will consider the case in which the group G is a Lie 
group which acts on H by a (projective) unitary repre- 
sentation {Ug}, whereas the initial state — also called seed 
state — is a pure state |^). Notice that the correspon- 
dence between transformed states and group elements is 
generally not injective, since the state j^E') may have a 
nontrivial stability group, say Gqr (we say that a group 
element h belongs to the stability group of j^*) when 
[//il^') = e^^'^l^*), with iph a real phase). In this way 
the transformed states are in one-to-one correspondence 
with the cosets gG^s, : in other words the group-orbit man- 
ifold (obviously invariant under the group representation 
{Ug}) is identified with the coset space X = G/G*. We 
see that in principle from the output state J7g|^) it is pos- 
sible to estimate the group element g of the transforma- 
tion Ug only if the stability group Gif of the input state 
I^E*) is trivial. Otherwise, we can estimate the coset x £ X 
which is in one-to-one correspondence with the output 
state I'^x) = 9{x) labeling any element of G in 

the coset x. In the following we will denote by xo = eG^ 
the coset containing the identity element e, and the seed 
state is relabeled accordingly as l^^o) = 1^)- This no- 
tation makes explicit the isomorphism between the coset 
space X and the homogeneous manifold of states \'^x) 
X £ X, i.e. on which the group acts transitively through 
its unitary representation as Ugl'^x) oc \'^gx) (apart from 
a phase factor). In this way, the estimation of the param- 
eter X € X becomes equivalent to a problem of covariant 
state estimation, and it was proved [ pO[ that the optimal 
probability distribution p{x\xo) of estimating x for in- 
put state I'i'xo) satisfies the identity p{gx\gxQ) — p(x\xq), 
namely the probability distribution on the manifold X for 
an input state Ug\'^) is equal to the probability distribu- 
tion for input state \'^) but with the manifold shifted by 
g^^. In the following we will suppose for simplicity that 
the group G is unimodular (i.e. the left invariant mea- 
sure d g on G is also right-invariant) and the stability 
subgroup is compact. According to a theorem by Holevo 
pO| , for square- integrable representations the covariant 
estimation is described by a POVM M on the probabil- 
ity space X of the general form 



dM(x) = C/,(.)SC/](^)d 



(1) 



where d x denotes the invariant measure on X induced by 
invariant measure d g on G p7| , and the positive kernel 
operator S belongs to the commutant G^ of the stability 
group (i.e. [S, = 0] , V/i G G*), and satisfies the 
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completeness constraint 

/ da;C/,(.)SC/;(^) ^ / dgU.EUl = /. (2) 

The fact that S e guarantees that the POVM does 
not depend on the particular choice of g{x). 



III. GROUP INTEGRALS OF OPERATORS 

The completeness constraint in (||) becomes particu- 
larly simple with some abstract considerations on group 
integrals. Since the group G is unimodular, its unitary 
square-summable representations satisfy Schur's lemma 
for any (generally infinite dimensional) representation 
space H Q, namely: 

For any couple {t/g } and {Ug} of irreducible com- 
ponents of the representation with invariant sub- 
spaces H^, H,y C H, respectively, every operator 
0^1/ : ^ satisfying the identity Uj^ O^i, = 
O^pUg Vg G G must be of the form 



0, otherwise, 



where ~ denotes equivalence of irreducible compo- 
nents, fc is a constant, and /^i, : ^ is the 
isomorphism mapping the two equivalent compo- 
nents, namely L/^ = IiJ.iyUglli, V.g E G (/^^ is the 
orthogonal projector onto the invariant irreducible 
subspace H^). 

A simple consequence of Schur's lemma is the Wedder- 
burn decomposition of operators O such that Tr[/pyO] < 
oo V/i, v [M 



(3) 



Taking the expectation values of both sides of Eq. (j^) 
on an arbitrary element |el^^) of an orthonormal basis 
{|em^)} for one has 



d.gTr[(t/t|e(f))(e(f)|;7,)0], Vn . (4) 



Applying now the Wedderburn decomposition to the 
group average of projectors lei'''') (el'''' | and using invari- 
ance of the subspace H^, one obtains 



G 



(5) 



where is a constant to be evaluated. We then have 



(6) 



where can be determined by taking the expectation 
value of both sides of Eq. ^ on any normalized vector 

in H^, in particular on the vector jei'^''), leading to 



(7) 



G 



On the other hand, if the representations /i and v are 
equivalent, there are two orthonormal basis {lelf"*)} and 
{|em'')} for and Hi,, respectively, such that I^^, = 
J2n |en''')(en''''|. Now, taking the matrix element of both 



sides of Eq. (H) between vectors |en^^) and \en') one has 



d.gTr[C/t|eM)(e(^)|{/,0], (8) 



The invariance of both subspaces and H,^ along with 
Schur's lemma give the identity 



(9) 



for suitable constant b^^ to be determined. Upon substi- 
tuting the last equation into Eq. (ph gives 



(10) 



and the constant b^f^ can be determined by taking the 
matrix element of Eq. (^ between vectors |e[,'^'') and 
jei''''), namely 



b,,= / d5(ei'')|[/t|e(r))(e(r)|t/,|e('^)) 



(11) 



G 



Notice that for equivalent components ^ ^ v for our 
choice of bases one has {e'n\Ug\e"n'') = {en''\Ug\en'') ^ 
whence b^^ 
the decomposition 



b^fj, = b^ij = b^. Summarizing, we have 



I dg UgOUl ^Y.^^J2 Tr[/,,,0]/^, , 
dg\{ei^^\Ug\el^^)\\ 

G 

If the group G is compact and its measure d g is nor- 
malized (i.e. J(^dg = 1), then it is easy to show that 
6p = where = dim(H^) (irreducible represen- 
tations of compact groups are finite-dimensional). In 
fact, summation over all n in Eqs. (^) and (|^) pro- 
vide in a direct way the values a^^ = Tr[I^^O]/dfi and 
a^^ = Tr[J^yO]/d^ for the coefficients in Eq. (§). On 
the other hand, the derivation given above holds for uni- 
tary square-summable representations, even with Dirac- 
orthogonal basis {|ei'^^)} for H^, namely (e^|e^,) — S{x — 
x'). The coefficients b~^ are generally non-integer, are 
called formal dimensions, and carry information about 
the structure of the irreducible components of the group 
representation. 
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IV. MEASUREMENTS WITH MAXIMUM 
LIKELIHOOD 

We will now consider measurements which maximize 
the likelihood, namely the conditional probability density 
p{x\x) of having the outcome equal to the true value for 
any x. Because of covariance this optimality criterion is 
equivalent to maximize the likelihood functional Cp [S] — 
Tr[2p] with p = |*)(vE'|, l^f) being the input state. 

Notice that the general solution to the maximum likeli- 
hood problem, which at first sight may appear of limited 
value, is actually equivalent to the solution of any quan- 
tum estimation problem with positive summable "goal"- 
function f{x,x) (the "goal" -function is the opposite of 
the customary "cost" -function —f{x,x) 0). This con- 
sists in associating to each measurement outcome x a 
"score" f(x,x), with the function f{x,x) increasing ver- 
sus X for X approaching the true value x. Then, the op- 
timal measurement is the one which maximizes the aver- 
age score. In a covariant estimation problem a mean- 
ingful goal function must satisfy the invariance prop- 
erty f{x,x) — f{gx,gx) \fg g G, and this allows to 
define a function h(g,g) on the group via the relation 
h{g,g) = f{gxo, gxo) for fixed xq. Then, the function h is 
positive (bounded from below), summable, and satisfies 
^{9t9) — h{g~^g,e), e denoting the identity transforma- 
tion. Now, thanks to covariance the average score can be 
written as 



dg%,e)Tr[p[/gSC/t] 



G 



dgh{g,e) 



C 



where 



^(p) 



J^dgh{g,e)UlpUg 
J^dgh{g,e) 



is a completely positive trace preserving map. Therefore, 
the maximization of a goal function can be viewed as 
a maximum likelihood scheme on the transformed state 
^{p), and depending on the form of the function h the 
choice of the input state may be restricted to special 
states, possibly mixed. Nevertheless, in this paper we 
will give a complete solution only for pure input states. 

The problem is now to find a positive operator S which 
maximizes the likelihood functional £p[S] = Tr[Sp], and, 
at the same time, satisfies the completeness constraints 
(^. Once an optimal S is found, the presence of a non- 
trivial stability group for l^t) can be taken into ac- 
count by replacing S with its group average over G>p 



_^ !^^dgU,~Ul 



(13) 



Notice that the value of the likelihood functional remains 
unchanged after this replacement, and the group average 
is still optimal (it is easy to show that the same occurs 



with ^{p) in the case of a general goal function). As a 
consequence of the Wedderburn decomposition (|l^), the 
completeness constraint (H) for S can be written as 



Tr[VS] 



(14) 



It is now convenient to decompose the input state \^) 
over the invariant subspaces of the representation as 
l^*) — X^^c^I^m)- This allows us to simply derive the 
following chain of inequalities 

= ^c;c,(*^|s|vi',) |c^||c.||Cm.I 

where the sums range in the set Mij, of all invariant sub- 
spaces which are nonorthogonal to j^"), denotes 
the likelihood functional defined by the pure state j^*), 
and denotes the matrix element (^p}^\^u). The 
first inequality can be saturated by the choice = 
.^yj^gj-g jg phase of c^. The second 
inequality is a necessary condition for positivity of S, 
and saturates for = yj (,^ni^vv (notice that this in- 
equality is not also a sufficient condition for positivity, 
whence the positivity of the optimal S must be checked 
a posteriori) . The third inequality is due to the fact that 
^^iM ^ Tr[/p^S] = h'p^ . Finally, the last Schwartz in- 
equality sets the following general upper bound for the 
maximum likelihood of covariant measurements 



-1 



(15) 



In the case of a compact group the inequality ( [T^ ) implies 
that the likelihood is always less than the sum of dimen- 
sions of invariant subspaces supporting \^). For infinite 
dimensions, on the other hand, the bound (|l5|) and the 
likelihood itself may diverge. One can see now that the 
following choice of the operator S 



(16) 



attains the bound (X^^eM* Ic/ily ^m^)^ fo'^ the likehhood 
functional. Note that, if \^) has no component in some 
irreducible subspace H,^, then the operator S must be 
extended to the whole space H, in order to fulfill the 
constraints Tr[/^^S] — for all /i. Obviously, such 
extension is generally not unique, e.g. one can take 



(17) 



where is any normalized vector in H,y, which both 
guarantees S > and satisfies the constraints Tr[/ppS] — 
for all /i. Notice that the presence of equivalent rep- 
resentations in Eq. generally improves the likelihood 
(this feature was missed in Refs. 16, EtI fsl). 



5 



If there are no equivalent representations in the de- 
composition of l^*), then the kernel averaged over 
the stability subgroup of Ivf) is optimal. However, in 
the presence of equivalent representations, one also wants 
the off-diagonal constraints Tr[/p,yS] = to be satisfied 
V/i ^ v. One can see that the kernel in Eq. satisfies 
also the off-diagonal constraints when the decomposition 



I*) 



satisfies 



(18) 



As shown in the Appendix, the subspaces carrying equiv- 
alent irreducible components of the representation can 
always be chosen in such a way to satisfy Eq. (18). It 
is worth noticing that the present "canonical" form for 
maximum likelihood measurements generalizes the case 
of the optimal covariant phase estimation given by Holevo 
pO| , further generalized in Ref. Finally, notice that 
the result derived here also holds for discrete groups, such 
as the permutation group or x by just substituting 
integrals with sums. 



subspaces than just those belonging to the symmetric 
space. Finally, the breach of the semi-classical limit in 
Eq. (^9|) does not necessarily need entangled POVM el- 
ements, and, actually, in one of the following examples 
the linear limit is overcome by a separable POVM. 



A. Universal state estimation 

1. SU{d) -covariant estimation: pure state estimation 

The estimation of a pure state in a finite dimensional 
Hilbert space H can be regarded as a covariant estimation 
with respect to the defining representation of the group 
SU{d), where d = dim(H). Indeed, the orbit of a given 
pure state contains all pure states of H. Clearly, the 
optimal kernel is S = d\4>){ip\, according to Refs. [g[ po[ , 
and consequently the value of the maximum likelihood is 
£(1) = d. 



EXAMPLES 



Pure state estimation with N > I copies in the same 
state 



While it is obvious that averaging the result over a 
number > 1 of equally prepared identical copies always 
improves the precision of estimation — either classically or 
not — a legitimate question is whether non-independent 
measurements on copies can be exploited to further en- 
hance the precision, compared to this conventional in- 
dependent measurement scheme. For the maximum- 
likelihood strategy, when measurements are performed 
independently on each copy, the value of the likelihood 
is linearly bounded as follows 

= /.../ dxi • • -dxNPixilx) ■ ■ -pixwlx) 

J J (19) 

-(^-^) 

< Nmax{p{xN\x)} = NC^'^^ , 

where is the maximum likelihood of a single-site 
measurement. 

As we will show in the following, the optimal measure- 
ment on N copies of the same state can surpass the linear 
bound for this semi- classical scheme involving indepen- 
dent measurements. Moreover, if we relax the restric- 
tion of identical preparation, corresponding to an input 
state of the form j'l') = jf/;)®^, it is also likely that a 
preparation in different states can lead to further im- 
provement in the estimation of the group-transformation 
Ug that occurred on the input state, since the decompo- 
sition of the global state |^') may involve more invariant 



This corresponds to the case of estimation of the group 
element g € SU{d) in the reducible representation U®^ 
with initial state \'^) = \tp)'^'^ . There are inequivalent 
components corresponding to the symmetric subspace 
(H®^)^, along with all other permutation invariant sub- 
spaces. Since j^*) belongs to the symmetric subspace 
(H®^)+ only, the optimal S is not unique, e.g. we can 
take S = dsi\i^)W)^^ + Iw, where ds = C^^"^) is 
the dimension of (H®^)+ and W is the orthogonal com- 
plement of (H'^''^)_|_. In any case with A^ copies we have 
maximum likehhood = dim(H®^)+, and for A > 2 
the semi-classical limit Nd is breached. Notice that the 
POVM is not separable, due to the presence of the or- 
thogonal projector /w. 



3. SU(d) estimation with two copies in different states 

In this case j^*) = \4')\<l>) can be decomposed as 
+ yin^lM/.), where s = KV'I^)! and 
l*±) = /J+ ± m^))- Then the optimal 

kernel ^ is proportional to the projector onto the vec- 
tor = y37|*+) + yrfTl*-) and the likelihood takes 

the value (y^ '^+'^+^'> + ^ '''^V^ )' ^ (by Schwartz 
inequality). It is easily seen that this bound can be 
attained by choosing — 2- The optimal POVM is 
separable (the optimal kernel is actually factorized as 
S — (i^ !?/')('!/' I ^ No further improvement can 

be achieved, since a likelihood greater than d^ is not 
compatible with the completeness of the POVM (in fact 
U[E]^TT[E] = d^). 
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B. Weyl-Heisenberg covariant estimation 

1. Estimation of displacement on the phase space 

This case corresponds to consider the Weyl-Heisenberg 
irreducible representation {D(z)} of the translation 
group on the complex plane, D{z) denoting the displace- 
ment operator D{z) = ~^ ° with [a, a^] = 1. Being 
non-compact, the representation space H is infinite di- 
mensional. Physically D{z) represents a joint shift of 
position and momentum of a quantum harmonic oscil- 
lator, and the covariant state estimation corresponds to 
a joint measurement of position and momentum. Here 
one has b — J^^-^\{n\D{z)\n)\'^, where \n) denotes an 
element of any orthonormal basis for H, which we can 
conveniently take as the set of eigenstates of the number 
operator a^a. Choosing n = one obtains 6=1, whence 
the optimal kernel for initial state {ip) is E = {4>\ and 
the maximum likelihood is £[S] = 1. Notice that for 
= |0) we get the well-known coherent-state POVM 
describing the heterodyne measurement O, EJ] . 



2. Estimation of displacement with identical shifts on 
N > 1 quantum oscillators 

This case corresponds to the tensor representation 
{D^^{z)} of the Weyl-Heisenberg group. The irre- 
ducible representations can be easily obtained by the lin- 
ear change of modes represented by the unitary transfor- 
mation 



where C„ 



V 



_ p'i>N[a[{a2-i halv)-al(a2^ ^'^jv)l 



with (h = 



■T^arctanx/iV^ so that VD®^{z)V'' = 
D{\/Nz) ® j®{N-i) ^ Then the irreducible subspaces are 
given by H„ = {V^W) ® \^n), W) & H}, where {|$„)} is 
an orthonormal basis for H*^'^"^^. The formal dimension 
coefficients are easily obtained as follows 



d z , 



IT 
j2 



|(o|($„|ra®^(z)yt|o)|ci>„)|2 

/ ^|(0|($„|i?(ViVz) 



1$. 



1 

TV 



Since the invariant subspaces carry all equivalent rep- 
resentations — the isomorphism between two of them is 
Imn = V^^(/(8) \'^m){'^n\)V — the problem of choosing 
a suitable decomposition of the initial state \^) in ir- 
reducible representations arises. In the general case, 
one should apply the full construction showed in Ap- 
pendix, while a simpler solution is possible for states of 
the form \^) = |ji)|i2) • • ■ |«jv)- In this case, one has 
only to write V\^) = E»i,»2,...,»„ c^l^2...^« l«i> 1*2) •• ■ Vn), 
and to define |<I>„) = X]j2...ijv '^™2...iN N2) • • ■ Nat) 



obtaining the de- 



sired decomposition j^I') = En CnV^\n)\<^n) (notice that 
{'^m\'^n) = Smn sincc they are eigenstates corresponding 
to different eigenvalues of the number operator). The 
value of the likelihood is then C[E] = A^(En ^n)^- 
We now consider two special cases. 

i) N copies of vacuum state |0): this case corresponds to 
the estimation of the complex shift z on the set {l^;)®^} 
of N copies of a coherent state \z). Here, the vacuum 
state l^") = |0)®^ belongs just to one invariant subspace, 
since Vl"^) — j^"). The optimal kernel is not unique, and 
IS given by any completion of S = iV(|0)(0|)^. The like- 
lihood value is N, namely the " semi-classical" value. For 

= 2 it can be proved that one of the optimal POVM's 
corresponds to averaging the outcomes of independent 
heterodyne measurements on two copies, while another 
optimal one corresponds to the independent measure- 
ment of the position ^{ai + o,\) and the momentum 

^{0,2 — 0-2)' taking as the outcome a = x + iy, where 
X and y are the two separate outcomes. 

ii) Two copies of a number state: |^) — \n)\n) with 
71 > 0. The maximum value of the likelihood is >C[S] = 



2 (ELo 2^(fc)v/(2fc)!(2«- 2fc)!j , and numerical cal- 
culation up to n = 1000 shows an almost linear be- 
versus n, much better than the semi-classical 
In the case of two copies of a one-photon state 
|1)|1). Decomposing the seed state we obtain 
-^(^■^120) +V^02)): an example of optimal 



havior 
value. 

I*> = 



kernel is then S = 2 2 (|1)(1|)®2 + E,-,o,2 V^0i){0i\V 
achieving a likelihood equal to 4, which is yet twice the 
semi-classical value. 



3. Estimation of displacement on two copies, with identical 
shifts in position and opposite shifts in momentum 

This case corresponds to the representation = 
D{z) ® D{z*)}, which is reducible, but does not pos- 
sess any irreducible proper component in H®^, and thus 
is beyond the hypotheses of our general results. In 
fact, the irreducible representations are all inequiva- 
lent, and make a continuous set, each component be- 
ing supported by the Dirac-normalized eigenvectors pl| ] 
-^\Diw)) - -^J2rn,niMD{w)\n) \m)\n) of the nor- 

mal operator W = I — I (^a^ (the heterodyne photo- 
current jl^, |lj, |32j). Upon expanding the operators 
V(z) = exp(zW^ — z*W) over the Dirac-orthonormal ba- 
sis, one has 



j2 

d z 



Viz)OV^iz) = / 



j2 

d z 



j2 

d w 



j2 / 

d w 



\Diw)){Diw)\0\Diw)){D{w)\ , 



j2 

d w , 
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namely a continuous version of the Wedderburn decom- 
position still holds 

f ^V{z)OVHz)= f d'wa^P^, 
Jc Jc 

for any O such that Ti-[P^O] < oo, with — 
\D{w)){D{w)\ and — n^^ Tt[PwO] (in a proper math- 
ematical setting the integral over w in the last equation 
should be interpreted as a direct integral). The maxi- 
mum likelihood covariant measurement for state estima- 
tion among the set generated by the seed j^") € H*^^ is 
given by the entangled kernel S = |??)(??|, where 

Jc I" 

which is the analogue of Eq. ( |l6|) for continuous spectrum 
(as in that previous case, 9^ is the phase of {D{w)\'i')). 
It is worth noticing that for — |0)|0) the problem 
corresponds to estimating the amplitude z of the set of 
coherent states and the value of the likelihood 

for the optimal measurement is 4, namely twice the like- 
lihood for the amplitude estimation for identical states 
{|z)|z)}. The probability distributions are indeed Gaus- 
sian in both cases, but the variance in this case is half 
the variance of the Gaussian for the states Iz)!^). The 
fidelity of the estimate is 2/3 for the states while 
it is 4/5 for |z)|2*). The last example can be regarded as 
the "continuous- variables" analogue of the measurement 
of the direction of two antiparallel spins by Gisin and 
Popescu [ p6[ . 

4-. Estimation of displacement on one part of a bipartite 
entangled system 

We consider here the representation {D(z) (g) /} acting 
on two optical modes. In this case the invariant sub- 
spaces are H„ = (g) \ipn),\ip) € H}, where {\(pn)} 
is any orthonormal basis in H, all of them support- 
ing equivalent representations, and the formal dimen- 
sions bn are all equal to 1. If we take a twin beam 
l^*) = \/l — x'^ X^ra 2;"|n)|n) as initial state, its decom- 
position is trivial, and are precisely the compo- 
nents l^'n) on the irreducible subspaces. Then the opti- 
mal POVM is given by \t]) = Y]^n)jn)^ame\y it is the 
two mode heterodyne POVM |13[ |3l| , p^ . Correspond- 
ingly, the value of the likelihood is £[S] = showing 
a strong enhancement by the effect of entanglement in 
agreement with Ref. 

VI. CONCLUSIONS 

By group theoretic arguments we have derived the class 
of measurements of covariant parameters which are opti- 
mal according to the maximum likelihood criterion. The 
optimization problem has been completely resolved for 



pure states under the simple hypotheses of unimodularity 
of the group and measurable stability group. The general 
method has been applied to the case of finite dimensional 
quantum state estimation with many input copies, and, 
for infinite dimensions, to the Weyl-Heisenberg covari- 
ant estimation, also giving a continuous-variables ana- 
logue of the estimation of direction on two antiparallel 
spins by Gisin and Popescu. The analysis allowed us 
to compare the "semi-classical" statistical scheme of re- 
peated identical measurements on identically prepared 
copies with non-independent measurement schemes, such 
as the LOCC and the entangled schemes. We have seen 
that for the maximum likelihood criterion, the possibil- 
ity of surpassing the semi-classical statistical efficiency 
with the number of copies N is essentially related to two 
factors: i) the way in which the dimension of A^-fold ten- 
sor product Hilbert space increases versus N; ii) its de- 
composition in irreducible subspaces. Moreover, the non 
separability of the optimal measurement — either in its 
POVM or in the optimal states — is strictly related to the 
structure of the group representation. 

We conclude by mentioning that the optimal covari- 
ant estimation for mixed input states is still an open 
problem, and an explicit analytical optimization seems 
a very difficult task. For the case of phase estima- 
tion the problem can be analytically solved in the spe- 
cial instance of states which are phase-pure For a 
general covariance group-representation the concept of 
phase-pure state can be generalized by choosing a vec- 
tor 1^';^) for each invariant subspace such that for ev- 
ery ^ ^ V one has = 0. Then every state 
p satisfying Supp{p} C Spanjl^-^)} and {^f,\p\^^) = 
g2(a;^-a:„)|^yj;^|^|^^^| bgi^^ygg ^s a purc statc in all re- 
spects, so that the upper bound and the canonical form 
of the optimal POVM still hold. It is likely that a gener- 
alization of this approach may extend the validity of the 
present solution of the covariant estimation problem for 
a special class of mixed states. 



VII. APPENDIX 

In this Appendix we show how to chose invariant sub- 
spaces of equivalent irreducible representations in or- 
der to satisfy Eq. (|l8|), namely = 0, for 
the decomposition j^E') — X^'^mI^m)' namely in such 
a way that all invariant subspaces effectively behave as 
supporting inequivalent irreducible representations. This 
choice of invariant subspaces guarantees that every op- 
erator S — 1 77) (77 1, where [77) is any linear combination 
of l^*!/), will satisfy the constraints Tr[S/^i,] — for any 
p ^ u. This method will allow to extend the general 
treatment of the phase estimation problem given in Ref. 
pof to any square-summable group representation. 

Let's consider an irreducible component of a unitary 
representation of G with multiplicity m < 00^ and de- 
note by H^") the invariant subspace carrying all equiva- 
lent irreducible components, and by H'^'^-' — a 
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given choice of invariant orthogonal subspaces H^, each 
carrying an equivalent irreducible representation. Since 
inequivalent irreducible components already satisfy Eq. 
(|l8|), wc can just focus attention on the component j^P^) 
of 1^-) on H^'^). 

Let's denote by /^^ the isomorphisms mapping into 
H^, and satisfying [I^u,Ug\ = Vg € G. As already 
mentioned, it is always possible to choose an orthonormal 
basis Bp — {lelf"*)} for each subspace in such a way 
that Bp = I^y^iy for any ^,v, where equality between 

bases is defined element- wise, i.e. \e[^'') — Ifiu\en'^) for 
all n. We have now the following simple lemma. 

Lemma (Choice of the decomposition into equivalent 
components). For each unitary matrix {V^i/} € the 
linear combinations B^ = X^i/^'i^Biy provide a new de- 
composition H'") = ®™^]^ Hp of H'") into subspaces sup- 
porting equivalent irreducible components, where = 
Span(B;,). 

Proof. The subspaces are orthogonal. In fact, 
upon defining B^ — {|//'^^)}, we obtain 

a, 13 

a 

Moreover, each carries a representation equivalent to 
that of, say Hi. In fact, the operator Sn = Ylv^t^'^^'^^ 
is indeed an isomorphism between the subspaces H i and 
Hp, since it defines a one-to-one correspondence between 
them via Bp = 5piBi, and commutes with Ug for all 
(7 G G, since each /pj^ commutes. This proves that 
the spaces {Hp} provide a new orthogonal decomposition 
H(") = ©pHp into invariant subspaces carrying equiva- 
lent components of the representation. 

Now, let's consider the component Y^S) of Y^) on H^'^\ 
and write its decomposition using the set of bases {Bp} 
as follows 

i*c.> = E*;;j4^V (20) 

We want to construct a new decomposition H'^, of H^") 
such that the components of on invariant subspaces 
satisfy Eq. (|l8|), namely they behave as belonging to 
inequivalent representations. This can be done as follows. 
Define oo > d = dim(Hp) and consider the m x d matrix 
^ {*J^^}. According to Eq. 5-" is a Hilbcrt- 

Schmidt operator, and hence we can write the singular 



value decomposition 

^'" = -l/^EC/, (21) 

where E is a m x d matrix with all vanishing off- 
diagonal elements, and V and U are m x m and d x d 
unitaries, respectively. From Eq. (|l]) one obtains 
— ^i^ti'^i^^i^iUin, where the sums run from 1 to 
r — rank(\I''^) < min(m,d). Equation ( pO[ ) then rewrites 

(22) 

where 

\9l''^)^J2v^,Ui^\ei'^). (23) 

fin 

The new bases B'^^ = {Iff;"^)} provide a new decompo- 
sition in invariant subspaces H'^ supporting equivalent 
components. In fact, starting from the set {Bp} the first 
unitary transformation U over each basis preserves the 
relations Bp = /p^B^, whereas the second transforma- 
tion according to the previous lemma, gives the new 
decomposition H^'^) = 0"'^^ H^, with H^ = Span(B;^). 
The state |*^) in Eq. ^ then satisfies 

(*plV|*j.) =0, for 1/ 

since the spaces H'^, = Span{|g[,'^'')} have been chosen such 
that each component of I'^^u) on HJ^ is just proportional 
to a single element l^^'^'') of the orthonormal basis with 
different I for different i^. 

Notice that, if the group is non-compact and there is an 
infinite number of equivalent irreducible subspaces, the 
spectrum of singular values of 4* may be continuous, and 
the sums in the above derivation must be replaced by in- 
tegrals, with some care in the generalization of definition 
and theorems. 
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